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Abstract 


T.R.  Bieler,  M.A.  Crimp,  D.E.  Mason* 
Michigan  State  University,  *  Albion  College 


An  integrated  experimental  and  theoretical  approach  to  measure  and  model  deformation 
transfer  and  microcrack  initiation  at  grain  boundaries  is  being  developed.  Electron  channeling 
contrast  imaging  (ECCI)  was  used  to  directly  image  deformation  defects  associated  with  strain 
transfer  and/or  crack  nucleation  at  grain  boundaries.  ECCI  observations  were  correlated  with 
mesoscale  crystallographic  texture  information  measured  by  selected  area  channeling  patterns 
(SACPs),  electron  backscattered  diffraction  (EBSD)  and  x-ray  texture  measurements.  To 
interpret  and  model  these  observed  grain  boundary  deformation  and  strain  accommodation 
phenomena,  a  nonlinear  variational  model  of  plastic  deformation  in  a  bicrystal  that  couples  bulk 
grain  deformation  with  compatibility  enforced  at  the  grain  boundary  was  partially  developed. 

Research  Objectives 

The  broad  objective  of  this  research  program  is  to  measure  and  model  the  process  of 
deformation  transfer  at  crystal  boundaries  in  anisotropic  materials,  with  a  view  towards 
predicting  optimal  microstructures  to  enhance  materials  for  Air  Force  applications.  A 
multidisciplinary  group  of  scientists  is:  1)  experimentally  measuring  and  quantitatively 
characterizing  dislocation  and  twin  structures  at  boundaries  in  polycrystalline  metals,  2) 
determining  how  crystal/boundary  misorientations  affect  deformation  transfer  across  boundaries, 
microcrack  initiation,  and  hence,  ductility,  and  3)  developing  quantitative  3 -dimensional  models 
that  describe  these  deformation  processes. 

Experimental  Approach 

The  experimental  research  objectives  were  pursued  using  polycrystalline  equiaxed  near 
gamma  TiAl  as  a  model  material.  Defect  generation  and  deformation  at  crystalline  boundaries 
was  examined  using  Electron  Channeling  Contrast  Imaging  (ECCI),  which  allows  imaging  of 
defect  structures  in  bulk  samples  using  an  FEG-SEM  and  channeling/difffaction  contrast 
[Simkin,  1999].  Two  loading  conditions  were  used.  Microcrack  nucleation  was  examined  at  a 
statistically  significant  number  of  grain  boundaries  under  a  known  global  state  of  stress  on  the 
tensile  surface  of  ex-situ  deformed  4-point  bend  specimens.  Crack  propagation  and  re- 
nucleation  was  studied  using  notched  4-point  bend  specimens  deformed  in-situ.  Observations 
were  correlated  with  the  grain-to-grain  misorientations  and  mesoscale  microtexture  determined 
using  electron  back  scattered  diffraction  (EBSD)  and  selected  area  channeling  patterns  (SACPs), 
allowing  determination  of  the  true  c  and  a  directions  in  the  tetragonal  TiAl  structure. 
Observations  of  cracked  and  intact  boundaries  were  analyzed  for  all  available  true  twinning  and 
ordinary  dislocation  systems  (neglecting  superdislocations)  in  grains  on  either  side  of  the 
boundary,  using  the  Schmid  factor,  the  Burgers  vectors,  and  the  stress  axis,  to  obtain  a  fracture 
initiation  parameter  that  can  be  used  in  computational  codes  to  predict  damage  nucleation. 
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Summary  of  Results 


This  is  the  final  report  for  work  conducted  from  12/14/00  to  12/13/03,  with  a  no  cost 
extension  to  8/13/04.  During  this  time,  two  Ph.D.  degrees  were  earned  by  B.  A.  Simkin  (2003) 
and  B.C.  Ng  (Ng  will  defend  his  dissertation  in  early  2005).  An  additional  graduate  student,  Y. 
Wang,  received  support  but  changed  projects.  In  addition,  visiting  scholars  K.  Boyapati  and 
Prof.  A.  Fallahi  were  supported  to  conduct  research  on  this  project,  and  an  M.S.  graduate 
student,  D.  Kumar,  was  supported  to  finish  analytical  work  that  followed  from  Simkin’s 
dissertation.  In  addition  to  the  $375,000  grant  from  AFOSR,  this  project  had  $1 14,621  in 
matching  funds  provided  by  the  MSU  Composite  Materials  and  Structures  Center. 

Brief  descriptions  of  major  accomplishments  are  provided  using  individual  researcher’s 
activities  as  a  focus.  Progress  reports  were  provided  in  July  of  2001, 2002,  and  2003,  and  copies 
of  oral  presentation  power  point  files  presented  in  August  contractors  meetings  have  been 
provided  earlier.  A  powerpoint  file  describing  progress  through  October  2004  was  submitted  on 
October  21,  2004. 

1.  Development  of  a  parameter  for  predicting  the  propensity  for  grain  boundary  microcrack 
nucleation  (B.A.  Simkin,  A.  Fallahi,  D.  Kumar) 

Using  the  experimental  paradigm  we  developed  in  this  program,  we  characterized 
deformation  mechanism  interactions  with  grain  boundaries.  In  particular,  we  found  that 
microcracks  nucleate  at  many  boundaries  where  deformation  twins  interact  with  grain 
boundaries  (Figure  1).  Close  examination  revealed  that  the  microcracks  can  be  attributed  to 
specific  twins  (see  dashed  arrows).  In  the  example  given  below,  there  are  two  twin  systems 
operating  in  the  upper  grain,  and  one  in  the  lower  grain. 


Figure  1:  ECCI  image  of  multiply  fractured  y-y  grain  boundary;  arrows  show  microcracks.  The  traces  of  {111} 
twin  planes  are  indicated  for  each  grain.  Ellipses  represent  unit  circles  on  slip  planes,  with  the  major  axis 
representing  the  plane  trace  evident  on  the  micrograph.  Arrows  represent  unit  vector  directions  for  ordinary 
dislocation  slip  (red)  or  twins  (dashed  black).  Bold  arrows  have  a  component  out  of  the  plane  of  the  page. 

Using  trace  analysis  based  on  SACP,  it  is  possible  to  determine  the  active  twinning 
planes  in  the  grains  on  either  side  of  a  boundary  (Figure  1),  because  the  orientation  of  the  grains 
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are  known  with  respect  to  the  macroscopic  tensile  stress  on  the  surface  of  the  ex-situ  deformed 
bend  samples.  Since  there  is  only  1  twinning  direction  on  each  twin  plane  in  TiAl,  it  is  then 
possible  to  completely  identify  the  twin  systems  involved  in  the  microcrack  nucleation.  The 
Schmid  factors  of  all  of  the  possible  twinning  and  dislocation  slip  systems  for  each  grain  were 
determined,  along  with  the  true  sense  of  shear  for  the  observed  deformation  twins  (the  twinning 
is  unidirectional).  The  slip  plane  inclination  is  shown  with  the  ellipses  (see  explanation  in 
caption). 

Though  the  Schmid  factors  are  moderately  high  for  all  three  of  the  twinning  systems 
identified,  the  microcracks  are  associated  with  the  twinning  system  on  the  (1-11)  plane  in  grain 
14,  which  has  the  lowest  Schmid  factor  of  the  three.  However,  this  is  the  only  twinning  system 
of  the  three  that  has  a  sense  of  shear  that  moves  mass  away  from  the  boundary  and  causes  an 
opening  Mode  I  type  stress  on  the  surface  at  the  grain  boundary.  This  tensile  mode  I  effect  on 
the  boundary  accounts  for  nearly  every  microcrack  studied.  In  every  other  boundary,  the 
microcracks  are  associated  with  the  twinning  system  that  has  the  highest  Schmid  factor  (Figure 
1,  where  the  twinning  system  associated  with  the  microcracks  it  is  the  third  highest  Schmid 
factor  is  the  only  exception).  Thus  microcrack  nucleation  is  most  likely  to  occur  on  boundaries 
that  have  highly  active  twinning  systems  that  move  mass  away  from  the  boundary  and  have  a 
high  Schmid  factor,  i.e.  a  high  probability  of  being  active,  and  hence  thicker  than  other  twins, 
making  the  sense  of  grain  boundary  strain  highest. 

In  a  recent  TEM  study,  Gibson  and  Forwood 
[2002]  found  that  twin  impingement  on  grain 
boundaried  in  y  TiAl  is  accommodated  by  ordinary 
dislocation  motion  on  both  sides  of  the  boundary. 

Based  on  this  observation,  Simkin  developed  a 
mechanistically  based  parameter  for  predicting  the 
propensity  for  microcrack  nucleation,  Fmax,  where  this 
factor  is  the  maximum  value  of  the  8  possible  values  of 
the  parameter  F  given  below  (one  for  each  of  the  4 
possible  twin  systems  on  each  side  of  the  grain 
boundary): 


F  =  mu 


b„„  t 


'tw 


St 

ord 


btw  "  bord 


Figure  2:  Populations  of  intact  and 
fractured  boundaries 


In  this  factor  mtw  is  the  Schmid  factor  for  a  specific  deformation  twinning  system  under  the 

/V 

global  stress  state,  the  bj  are  unit  vectors  in  the  sample  coordinate  system  that  describe  the 
Burgers’  vector  directions  for  the  each  twin  system  tw  and  all  8  of  the  ordinary  dislocation 
systems  in  each  pair  of  grains,  ord. 

The  rationale  behind  this  parameter  is  a  product  of  probabilities,  i.e.  F  identifies 
conditions  where  slip  transfer  is  highly  likely;  F  is  large  when  twinning  is  highly  likely  (large 

mM),  when  the  direction  of  b  M  has  a  significant  degree  of  alignment  with  the  tensile  axis  such 
that  a  mode  I  opening  force  is  generated  in  the  boundary,  and  when  the  twinning  vector  is 
reasonably  aligned  with  ordinary  dislocation  slip  directions  that  can  accommodate  the  twin  shear 

in  either  grain,  i.e.  cos  k  =  b^  -bord  as  illustrated  in  Figure  3.  This  parameter  was  evaluated  to 

find  the  maximum  value  of  Fmax  for  each  boundary  1 1  cracked  and  1 1  intact  boundarys  with 
microstructurally  similar  geometric  characteristics.  The  values  of  Fmax  for  the  the  two 
populations  were  examined  using  the  student  t-test,  showing  that  the  two  populations  are 
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statistically  distinct  with  probability  >  99.5%,  represented  schematically  in  Figure  2  (assuming 
normal  distributions).  The  fact  that  fractured  boundaries  are  correlated  with  a  high  degree  of  slip 
transfer  is  ascribed  to  leaving  residual  slip  vectors  in  the  boundary  plane  that  increases  boundary 
strain  energy. 

Since  F  consists  of  parameters  that  can  be  evaluated  in  plasticity  based  models  of 
microstructures,  this  result  suggests  that  Fmax  may  predict  the  propensity  of  a  boundary  to 
develop  microcracks  and  hence  provide  a  mechanistic  way  to  predict  damage  nucleation  in 
computation  models.  However,  in  the  form  given  above,  the  factor  does  not  take  in  to  account 
the  orientation  of  the  grain  bounday  relative  to  the  tensile  axis,  both  of  which  would  be  exected 
to  play  a  role  in  microcrack  nucleation.  This  improvement  can  be  expressed  by  components  that 
maximize  mode  I  opening  components,  where  n  is  the  grain  boundary  normal. 

^="Jtw|btw-“|t-n|E|^tw-bord 

ord 

This  parameter  was  evaluated  using  information  obtained  from 
serial  sectioning  to  obtain  grain  boundary  inclinations  by  Fallahi 
and  Kumar.  However,  knowing  the  grain  boundary  inclination 
has  did  not  improve  the  ability  of  F  to  separate  the  boundary 
populations,  which  was  surprising;  this  result  is  examined  in  more 
detail  in  a  paper  in  preparation  (Fallahi  et.  al). 

Deepak  Kumar  has  further  explored  a  large  number  of  _ .  ,  , 

variations  of  the  fracture  initiation  parameter  to  examine  what  transfer  at  grain 'boundary.*2111 

combination  of  factors  can  best  separate  the  cracked  and  intact 

grain  boundary  populations.  He  examined  how  misalignment  of  slip  planes  affects  the 
propensity  for  slip  transfer  which  was  not  considered  in  prior  models.  Ashmawi  and  Zikry  used 
a  parameter  to  quantify  this  misalignment,  that  can  be  expressed  as,  cos  9  =  n(r  -h(r  ord ,  the 

degree  that  the  traces  of  the  slip  planes  of  ordinary  dislocations  are  aligned  with  the  trace  of  the 
twin  plane  on  the  boundary  plane  which  has  been  used  as  and  additional  accomodation  factor  as 
follows,  and  this  version  of  F  has  been  the  most  successful  in  distinguishing  between  fractured 
and  intact  boundaries.  This  result  is  described  further  in  a  paper  by  Kumar  et  al.,  that  is  in 
preparation. 

■^2 max  =  maxK,  I  b*  t  1 1 ;  I  b  J  1 1 1  n,^  -n^  |) 


Elastic  anisotropy  also 
plays  an  important  role  on  the 
stress  in  a  grain  boundary 
region.  This  was  examined  in 
symbolic  triple  junctions  using 
the  known  three  dimensional 
stiffness  (C #*/)  in  each  grain,  as 
illustrated  in  Figure  4.  So  far, 

we  have  estimated  t  as  uniaxial 
tension,  but  the  traction  acting 
on  the  boundary  is  likely  to  be 
somewhat  different  from  the 
global  stress  state.  To  find  the 
local  stress  tensor  at  a  boundary 


Figure  4  Full  3-D  elastic  anisotropic  finite  element  analysis  of  von 


Mises  stress  distribution  in  a  TiAl  tri-crystal,  (calculation  provided  by 
A.  Zamiri  and  F.  Pourboghrat) 
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requires  elasto-plastic  modeling  of  actual  micro stmctures,  which  is  beyond  the  scope  of  this 
program.  With  knowledge  of  the  local  traction,  we  expect  the  distribution  depicted  in  Figure  2  to 
be  narrower. 

2.  Prediction  of  crack  path  using  a  fracture  propagation  parameter  (Ng) 

Ng  adapted  the  ideas  introduced  by  Simkin  to  develop  a  crack 
propagation  parameter.  He  identified  cracked  and  intact  boundaries  in 
the  notched  4-point  bend  specimen  in  Figure  5,  where  the  initial  crack 
propagated  and  arrested.  In  this  study,  the  role  of  twinning  as  a 
deformation  accommodation  mechanism  in  grain  B  due  to  twinning 
deformation  in  grain  A  was  examined,  because  crack  growth  involves  Aluminum 
strain  gradients.  This  requires  consideration  of  directionality,  since  slip  Figure  5  Notched  crack 

transfer  by  twinning  requires  that  the  absolute  direction  of  the  twinning  growth  specimen 
shear  is  not  opposed  to  twinning  in  the  neighboring  grain.  Directionality 

is  also  relevant  to  the  process  of  crack  propagation,  because  a  crack  approaching  grain  A  before 
grain  B  will  cause  more  deformation  in  grain  A  than  B.  Thus,  slip  transfer  is  more  important 
from  grain  A  to  B  than  B  to  A.  A  crack  propagation  parameter  can  be  expressed  as 

f  2  4 

^A->B  ~mAtw  \^Atw  ^  Bordy  Zj  ^  Atw  Btw 

\^Bord= 1  Btw- 1  if  >0 

The  three  parts  of  this  parameter  are  explained:  mAtw  is  the  maximum  Schmid  factor  for  a 
reference  twin  system  in  grain  A;  bAtw  t  identifies  how  well  the  reference  twin  Burgers  vector  is 

aligned  with  the  tensile  traction  unit  vector  t  (for  simplicity,  in  the  direction  of  tensile  loading); 
the  two  sums  of  the  dot  products  in  large  parentheses  describe  how  well  the  Burgers  vector  of 
this  most  highly  stressed  twin  system  in  grain  A  is  aligned  with  the  Burgers  vectors  for  ordinary 
dislocations  and  twins  in  grain  B.  Also,  using  a  clever  synthesis  of  EBSP  and  ECCI 
measurements,  Ng  was  able  to  characterize  crystal  orientations  for  more  than  twice  as  many 
grain  boundaries  as  Simkin.  The  histogram  in  Figure  6  shows  that  microcracking  is  unlikely 
when  Fa->b  <  1.0,  but  is  more  likely  when  FA->b  >  1.0.  From  a  t-test  analysis  (using  a  non-pooled 
and  unequal  variance  mean  comparison,  assuming  the  underlying  populations  are  approximately 
normal),  the  mean  value  of  FA->b 
for  the  cracked  population,  Pabc,  5 

exceeded  the  intact  population 
mean,  Pabi,  with  99.9%  j=  4 

confidence.  Thus,  this  parameter  | 

may  be  used  to  predict  a  path  for  %  3 

crack  propagation  (i.e.,  grain  o 

boundaries  with  Fa~>b  >  1.0  are  2  2 

<D 

prone  to  fracture  while  boundaries  f 

with  lower  values  are  not).  z  1 

Figure  7  shows  a 
comparison  of  the  mapped  0 

microstructure  ahead  of  the 
arrested  crack  with  the  specimen 

after  the  crack  was  grown  with  Figure  6  -  Histogram  of  intact  and  cracked  boundaries  using  FA.B. 


0.3  0.4  0.5  0.6  0.7  0.8  0.9  1  1.1  1.2  1.3  1.4  1.5 

Fracture  propagation  parameter 
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additional  loading.  Use  of  the  Fa->b  to  identify  strong  and  weak  boundaries  are  indicated  in 
Figure  7  with  black  lines  for  weak  boundaries.  The  effectiveness  of  this  parameter  is  illustrated 
by  the  fact  that  the  crack  path  connected  regions  of  weak  boundaries.  Also,  intergranular 
cracking  occurred  in  boundaries  that  did  not  have  a  Mode  I  opening  force,  but  were  weak.  It  is 
notable  that  strong  boundaries  prevented  crack  propagation  and  that  an  annealing  twin  boundary 
was  oriented  to  be  a  weak  boundary  where  the  crack  jogged  parallel  to  the  tensile  axis.  This 
shows  that  grain  boundary  character  can  also  be  defined  in  terms  of  its  slip  characteristics. 

It  is  remarkable  that  grain  boundary  deformation  character  on  the  surface  affects  the 
fracture  path  as  much  as  it  does;  suggesting  that  microcracking  on  the  free  surface  may  be  more 
likely  than  beneath  the  surface  (consistent  with  the  geometry  of  twins  that  move  mass  away  from 
the  grain  boundary  identified  above).  Some  3-D  Anisotropic  elastic  finite  element  simulations  of 
deformed  TiAl  tri-crystals  show  the  highest  stress  concentrations  at  boundaries  on  the  surface. 


\  l  *  * 

Figure  7  Crack  propagation  and  map  of  microstructure  ahead  of  initial  crack  tip.  Bold  lines  represent  weak 
boundaries,  and  dotted  line  represents  actual  crack  path.  Not  all  weak  boundaries  cracked,  see  Figure  6. 


3.  Modeling  of  Slip  Transfer  in  a  Bi-Crystal  (Mason) 

Darren  Mason’s  contributions  to  this  effort  have  been  important,  though  less  visible.  His 
time  available  for  this  effort  was  affected  by  his  moving  to  Albion  College  during  the  course  of 
this  project.  Much  of  the  analytical  work  needed  to  accomplish  the  tasks  described  above  were 
made  possible  by  Darren’s  development  of  software  that  accelerated  microscopic  analysis.  In 
addition  to  this  practical  contribution,  he  also  developed  a  foundation  for  modeling  slip  based 
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grain  boundary  deformation  in  the  form  of  a  variational  model.  No  further  progress  occurred  in 
this  effort  since  Summer  2003,  and  the  state  of  this  model  development  is  repeated  from  the 
2003  progress  report: 

To  begin,  we  summarize  briefly  the  kinematical/kinetic  structure  of  this  problem.  Let 
O®,  j  =  A,  B  correspond  to  two  bounded  regular  three-dimensional  regions  with  piecewise 
smooth  boundaries  that  represent  the  reference  configuration  of  each  grain  in  the  bi-crystal.  The 
deformation  of  each  grain  is  modeled  by  a  smooth  mapping  y®  so  that  the  domain  occupied  by 
the  deformed  grain  is  to®  =  y®(Q®).  The  material  is  endowed  with  a  Helmholtz  bulk  free 
energy  yr. = y/  (F®,  y®)  as  well  as  an  interfacial  energy  density  %(H,  p)  representing  an  energetic 
contribution  resulting  from  lattice  incoherence  /  distortion  at  the  grain  boundary.  Here  y®  and  p 
correspond  to  vectors  of  internal  variables  and  H=[F(^]'1  [F^] 
denotes  the  relative  deformation  gradient  evaluated  at  the  grain 
boundary  (cf.  [Cermelli,  Kinderlehrer]). 

Within  this  framework,  given  a  deformation  history  F ®(f) 
at  a  fixed  material  point,  we  require  that  the  defect  state  y  ®(f)  in 
the  material  evolve  along  a  path  of  minimal  work  W,  where 


Undeformed 
Bi-Crvstal  Halves 


W=  \^--F  dt+  ^-FLdt. 

J  dF  J  ^TT 


5H 


Integration  by  parts  at  any  material  point  yields 


d% 


dy 


so  that  calculation  of  8  W  =  0  at  a  point  away  from  the  grain 
boundary  for  variations  is  equivalent  to  that  found  in  [Ortiz]. 

On  the  other  hand,  consideration  of  points  on  the  grain 
boundary  requires  more  care.  One  approach  is  to  consider  the 
framework  outlined  in  Figure  8.  Here,  Bj(h)  represents  the  pre¬ 
image  of  a  neighborhood  of  the  boundary  point  in  reference  grain 
j.  The  energy  associated  with  that  point  then  becomes 


Figure  8.  Schematic  image  of 
reference  and  deformed  bi¬ 
crystal  configuration  for  grain 
boundary  model  framework. 
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Variational  analysis  of  this  energetic  expression  for  compatible  deformations  of  the  grains  yields 
the  equilibrium  equation 


SW  =  5WV  +  SWZ  =  0  = 
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where  G  ,  tj,  and  k  respectively  represent  arbitrary  admissible  variations  of  the  grain 
deformation  gradient,  bulk  and  grain  boundary  internal  variables. 
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The  experimental  determination  of  a  fracture  initiation  parameter  described  above 
suggests  that  boundaries  that  are  not  good  at  transferring  slip  are  less  likely  to  crack.  This 
physical  result  is  highly  compatible  with  the  variational  modeling  approach,  where  the 
deformation  necessary  to  maintain  the  boundary  compatibility  must  be  accomplished  by 
deformation  within  the  grain  and  not  across  the  grain  boundary.  When  this  variational  model  is 
embedded  within  a  homogeneous  surrounding  medium  deforming  under  the  influence  of  the 
boundary  conditions,  the  local  stresses  and  strains  in  the  boundary  region  near  each  grain  can  be 
evaluated  accurately.  This  model  will  generate  strains  near  the  boundary  that  will  differ  from  the 
grain  interior,  as  observed  in  microscopy.  Since  the  boundary  in  the  variational  model  has  an 
energy  associated  with  it,  it  will  be  possible  to  investigate  and  incorporate  fracture  energy 
criteria  that  lead  to  separation.  Thus,  this  variational  model  has  the  potential  to  provide  a 
physically  and  microstructurally  based  damage  initiation  condition  that  can  be  incorporated  into 
a  larger  scale  continuum  damage  mechanics  FEM  computation.  With  this  stage  of  development, 
practical  tools  for  design  of  materials  for  critical  applications  will  be  available  to  enhance  the  Air 
Force’s  ability  to  introduce  advanced  materials  with  confidence. 

Summary 

Our  investigation  of  microcrack  nucleation  has  resulted  in  identification  of  a  class  of 
fracture  initiation  parameters  that  are  based  upon  slip  activity  in  adjacent  grains.  This  parameter 
can  be  used  to  identify  conditions  for  damage  nucleation  in  computational  micromechanical 
models.  Thus,  we  have  identified  a  slip  based  definition  of  grain  boundary  character  that  is 
capable  of  predicting  damage  nucleation  and  the  crack  propagation  path  that  is  based  upon  the 
physics  of  micromechanical  deformation. 
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